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Abstract 

The  incremental  shear  modulus  is  calculated  for  infinitesimal 
shear  of  an  incompressible  isotropic  elastic  material  in  a state 
of  pure  homogeneous  deformation.  Necessary  and  sufficient  condi- 
tions that  it  be  positive  are  obtained.  These  conditions  are 
related  to  the  Hadamard  conditions  for  propagation  of  a wave  of 
infinitesimal  amplitude  in  the  deformed  material. 
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1 . Introduction 

In  the  present  paper  we  calculuate  the  incremental  shear 
modulus  for  infinitesimal  simple  shear  of  an  incompressible 
isotropic  elastic  material  which  is  initially  subjected  to  a 
pure  homogeneous  deformation,  the  direction  of  shear  and  plane 
of  shear  being  arbitrary.  (An  analogous  calculation  has  been 
carried  out  previously  [1]  for  the  case  when  the  plane  of  shear 
is  a principal  plane  for  the  pure  homogeneous  deformation.) 

Necessary  and  sufficient  conditions  are  obtained  for  the 
incremental  shear  modulus  to  be  positive  for  all  directions  of 
shear  perpendicular  to  a fixed  direction,  which  lies  in  the 
plane  of  shear  and  is  defined  by  a unit  vector  n say.  These 
are  the  inequalities  (4.8)  of  the  present  paper.  It  is  shown 
in  §5  that  they  are  equivalent  to  the  necessary  and  sufficient 
conditions,  previously  obtained  in  [2]  and  given  in  (6.5)  of 
the  present  paper,  that  the  two  sinusoidal  waves  of  infinitesimal 
amplitude  which  can  propagate  in  the  deformed  elastic  material, 
in  the  direction  n , have  real  velocities.  The  velocities  of 
propagation  of  these  waves  are  given  (see  equation  (5.6)  below) 
by  /y/p  , where  p is  the  material  density  and  y is  one  or 
other  of  the  stationary  values  of  the  incremental  shear  modulus 
for  shearing  in  a direction  perpendicular  to  n . It  is  shown 
in  §7  that  the  directions  of  polarization  for  the  waves  are  the 
directions  for  which  the  incremental  shear  modulus  takes  its 
stationary  values.  ! 

The  necessary  and  sufficient  conditions  (6.5)  for  the  wave 
velocities  to  be  real,  and  the  equivalent  conditions  (4.8)  for 

J 


the  incremental  shear  moduli  to  be  positive,  involve  the  first 


and  second  derivatives  of  the  strain-energy  function  for  the 
material  with  respect  to  the  strain  invariants,  as  well  as  the 
principal  extension  ratios  for  the  pure  homogeneous  deformation 
and  the  unit  vector  n . It  is  desirable  to  obtain  these  condi- 


tions in  a form  which  is  independent  of  n . This  has  been 
achieved  in  [1,3]  in  the  case  when  n is  restricted  to  lie  in 
a principal  plane  of  the  pure  homogeneous  deformation.  It  has 
already  been  shown  that  the  conditions  so  obtained  (the  inequali- 
ties (6.10)  of  the  present  paper)  imply  the  condition  ( 6 . 5 ) 2 . 

In  the  present  paper  it  is  shown  that  a further  one  of  the  condi- 
tions (condition  (6.5)^)  is  automatically  satisfied.  It  was  shown 
in  a previous  paper  [4],  by  numerical  example,  that  the  first  of 
the  conditions  (6.5)  is  not  implied  by  the  conditions  (6.10). 
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2 . Incremental  constitutive  equation 

We  consider  a deformation  of  an  incompressible  isotropic 
elastic  material  in  which  a particle  which  is  initially  at 
vector  position  £ with  respect  to  a fixed  origin  0 , say, 
moves  to  vector  position  x . The  components  of  £ and  x 
in  a fixed  rectangular  cartesian  system  x , with  origin  at 
0 , are  and  respectively.  The  deformation  gradient 

matrix  g is  defined  by* 


8 ■ IlgJI  * Hxi,all  • (2.1) 

The  Finger  strain  C is  defined  by 

£ - gg+  , (2.2) 

where  the  dagger  denotes  the  transpose.  Let  W be  the  strain- 
energy  per  unit  volume.  Then,  since  the  material  is  isotropic 
and  incompressible,  W must  be  expressible  as  a function  of 
the  invariants  T1  and  T2  of  C , defined  by 

Tj^-trC,  T2  « ^-{(trC)2  - tr  £2}  . (2.3) 

Thus, 


W - W(TlfT2)  . 


(2.4) 


The  Cauchy  stress  matrix  a * 
a * 2 1 


'U1 


is  then  given  by 

- :[(It1*ti  -2]'FS- ' (2-5) 

where  6 is  the  unit  matrix  and  F is  an  arbitrary  hydrostatic 
pressure. 


We  use  the  notation  „ to  denote  differentiation  with  respect 

> o 

to  £a  and  shall  assume  that  the  Einstein  summation  convention 
applies  to  lower  case  Greek  and  Latin  subscripts. 


» 
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We  now  consider  that  the  deformation  £ -*■  x is  the  resultant 
of  a finite  deformation,  in  which  the  particle  initially  at  vector 
position  £ moves  to  vector  position  X , followed  by  an  infini- 
tesimal deformation  in  which  it  moves  to  vector  position  x . We 


shall  accordingly  write 


and  correspondingly 


x ■ X ♦ u 


(2.6) 


O ■ o+  s , Tx  - Ix  ♦ , T2  • I2*  i2  , 

C • C * c , P ■ P ♦ p , 


(2.7) 


where  o,  I1,  I2,  C and  P are  the  values  of  o,  T , I2>  C and  P 
associated  with  the  deformation  £ -►  X . 

v — 

With  this  notation,  it  follows  from  (2.5)  that 
s - 2[(W1*I1W2)c  - W2(Cc+cC)  ♦ ixW2C 

tW12*IlW22)ij)?'  (W12il*W22i2)92]'  t2'8) 


where 


w - 3W 

"i  »T 


w . 3W 

2 iT, 

u«0  * u«0 


. ", 


, etc.  (2.9) 


1 u-0 


It  follows  from  (2.1),  (2.2),  (2.3),  (2.6)  and  (2.7)2>3>u, 
with  the  neglect  of  terms  of  higher  degree  than  the  first  in  the 
spatial  derivatives  of  u , that 

ij  i » ci  J * oi  J * oi  i » oi 

(2.10) 

ll  " 2Xi,aui,a  ’ Wl’Vji' 


(2.10) 


6. 


3.  The  incremental  shear  modulus 

We  now  consider  that  the  deformation  £ -►  X is  a pure 
homogeneous  deformation  with  principal  extension  ratios  X , 
X2 , and  principal  directions  parallel  to  the  axes  of  the 


rectangular  cartesian  coordinate  system  x . Thus  , 


Xa£a  (A-  1,2,3)  . 


(3.1) 


We  suppose  that  a uniform  simple  shear  of  amount  k is 
superposed  on  this  pure  homogeneous  deformation,  the  direction 
of  shear  being  parallel  to  the  unit  vector  t , say,  and  the 
plane  of  shear  being  parallel  to  the  unit  vectors  l and  n , 
say,  where  n is  perpendicular  to  t . Then,  in  (2.6)  u is 
given  by 

u - <(n -XU  • (3.2) 

With  (3.1),  equations  (2.1),  (2.2)  and  (2.3)  yield 


'AB 


*ASAB  • 


^ ’ l X 


B 


I, 


5*i 


(3.3) 


where  is  the  Kronecker  delta  and  T denotes  summation 

AB  g 

over  B-  1,2,3  . The  incompressibility  of  the  material  implies 
that 


X1X2X3 


(3.4) 


Introducing  (3.1),  (3.2)  and  (3.3)1  2 into  (2.10),  we  obtain, 
with  £*n»0  and  (3.4), 


We  shall  assume  that  the  Einstein  summation  convention  is 
not  applicable  to  upper  case  Latin  subscripts. 
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'AB 


*^AnBXB*2B  "a'a*  ’ 


li  - 2k  i Vax;  • 

A 

12  " XB  ]•  ^AnAXA  " | ^AnAXA^  2lC  | ^AnAX A^  * 

From  (3.3).^  and  (3.5)  , we  obtain 

(?2)ab  - *i«AB . 

<S9+CC)AB  - <(Va\2+Vb^(^^  • 


(3.5) 


(3.6) 


Introducing  (3.3)1  2,  (3.5)  and  (3.6) 

SAB  " 2k*AB  * P<SAB  * 


into  (2.8)  we  obtain 
(3.7) 


where 


?AB  ■ WBXB*  VaXA> 

* VVb\2*  Vaxb2‘“ab  I Vc‘^1 

* 2xa5abCku  l Vcx?‘"i2<(VXA>  l Wc 

■ | Vcxc2l  • I Vcxc2]  • 


(3.8) 


In  deriving  this  result,  we  use  (3.4)  and  the  relations 

i t-l  - i n!  ■ i . i v*  ■ 0 • (3-9) 

We  note  that  the  plane  material  surface  which  is  normal  to 
the  unit  vector  n after  the  pure  homogeneous  deformation  remains 
normal  to  n after  the  imposition  of  the  simple  shearing  deforma- 
tion and  the  area  of  an  element  of  this  surface  is  unchanged  by 
this  deformation.  Accordingly,  the  change  f in  the  component 
of  the  stress  vector  acting  on  this  surface  in  the  direction  of 
the  unit  vector  t is  given  by 


W 


8. 


“ SAB^AnB 


- SAB£AnB  ' 


A , B " “ A , B 

The  incremental  shear  modulus,  u say,  for  the  simple  shear 
considered  is  given  by 


(3.10) 


u ■ f/< 


(3.11) 


Then,  introducing  (3.8)  into  (3.10)  and  using  the  relations 


(3.9)  and  (3.5)^,  we  obtain 

T?  ’ "l  l XInA*W2  l XA2( 


2,2 

A 


* Wi)‘-!»u({  xaVami  x;2Va) 

* **22  0 XA  {AnA)  3 ' 

A 


(3.12) 


9. 


4 . Conditions  for  the  incremental  shear  modulus  to  be  positive 
In  this  section  we  consider  that  the  direction  of  the  unit 
vector  n is  fixed  and  we  derive  conditions  for  the  incremental 
shear  modulus  to  be  positive  for  all  choices  of  the  direction  of 
shear  l perpendicular  to  n . 

Let  a,  b be  two  perpendicular  unit  vectors  in  a plane 
normal  to  n . We  choose  their  positive  directions  so  that 
a,  b,  n forms  a right-handed  triad.  Then, 


a*a  ■ b*b  * n*n  ■ 1 , a*b  ■ a*n  » b*n  ■ 0 , n ■ a*b  . (4.1) 
Since  t is  perpendicular  to  n , we  may  express  it  in  the  form 


t ■ a cos  9 ♦ b sin  9 


(4.2) 


Introducing  (4.2)  into  (3.12),  we  obtain 


where 


jU  - P cos  29  ♦ Q sin  29 ♦ R , 


(4.3) 


P - W2P1^2W11P2-4W12P3+2W22PU  , 

Q - W2Q1+  2WllQa-  4W12Q3+  2W22Qu  , 

R - WiRq  ♦ W2Rx  ♦ 2WX1R2  - 4W12R3  ♦ 2W22RU  , 

and  PltP2, . . . .R^  are  defined  by 

P1  » j l *A2(aA"  bA5  ’ R1  • 7 I Xa2(®a+  bA>  * 
A A 

<5i  ' I *A2aAb»  : 

A 


(4.4) 
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— 


►. 

^2  " XAnAaA^J  XAnAbA^  ’ 


r2 

x*v/* 

7^  XAnA! 
A 

xa!2  • 

P3 

■ if  I 

4 A 

^VaHI 

A 

XA2nAaA) 

Xf"AbAHJ 

A 

XA2nAbAf 

XaVaHI 

A 

XA"AbAH^ 

A 

XA2"AaAf 

R3 

XAnAaA  f f ^ 

XI2nAaA> 

**{I 

XA°AbA  f f ^ 
A 

x;VaI 

pu 

' ^ A*2' 

> 

■ (I  * 

A 

A 

XA2"AbAl 

) 

- p * 2 

XA  nAaAt 

* x;2' 

,2 

’AbA> 

» 

Ro 

* I XA 

A M 

"a  • 

We  can  rewrite  (4.3)  as 

ju  “ (P2+Q2)issin(26+x)  ♦ R , (4.6) 

where 

X - tan"1 (P/Q)  . (4.7) 

It  is  apparent  that  the  necessary  and  sufficient  conditions  for 

y to  be  positive  for  all  values  of  8 are 

R > 0 and  R2  > P2  + Q2  . (4.8) 

The  stationary  values  of  u occur  when 

9-^-|  (m-  1,3, 5, 7)  . (4.9) 

m»  5,7  correspond  to  reversal  of  the  directions  of  L corre- 
sponding to  m*  1,3  respectively.  We  need  therefore  consider 


11 . 


i 

only  the  cases  when  m=  1,3  . The  two  corresponding  directions 
of  i are  perpendicular.  Whem  ni'l  , u has  a maximum  value 
and  when  m-  3 , u has  a minimum  value.  From  (4.6)  it  is 
evident  that  these  values  are  given  by 

U - 2{R±  (P2+Q2) . (4.10) 


I 


12. 

5 . Propagation  of  a plane  sinusoidal  wave 

We  now  suppose  that  instead  of  subjecting  the  pure  homoge- 
nously deformed  material  to  a static  infinitesimal  simple  shearing 
deformation,  we  propagate,  in  the  direction  n , a cosinusoidal 
shear  wave  of  infinitesimal  amplitude  < and  angular  frequency 
oj  , polarized  in  the  direction  l . We  note  that  for  an  incom- 
pressible material,  the  directions  of  propagation  and  polarization 
are  necessarily  perpendicular. 

With  the  usual  complex  notation,  the  displacement  field  u 
associated  with  this  wave  is  given  by 

u = <1  exp  iu)(Sn*X-  t)  . (5.1) 

It  is  evident  from  (2.8)  and  (2.10)  that  the  incremental  stress 
associated  with  the  deformation  (5.1)  is  obtained  from  (3.7)  by 
replacing  < by  ;<jjSk  exp  iu)(Sn*X-  t)  , ‘thus 

sAB  = 2iu)Sk?ab  exp  ia)(Sn*X  - t)  - p6AB  ; (5.2) 

where  ?AB  is  given  by  (3.8).  We  take  the  incremental  hydro- 
static pressure  p in  the  form 

p = p exp  iai(Sn*X-  t)  , (5.3) 

where  p is  independent  of  X and  t . 

Bearing  in  mind  the  fact  that  the  deformation  £ •+■  X is 
homogeneous,  we  obtain  the  incremental  equation  of  motion  as 

| TV  ' c“a  ’ t5-4) 

where  p is  the  material  density. 


1 
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Introducing  (5.1),  (5.2),  (5.3)  into  (5.4),  we  obtain 


nBSAB 


♦ JJL  n - 4 r 

<i»)S  A g2  J 


Multiplying  (5.5)  throughout  by  t ^ and  summing  over 
we  have,  with  (3.9),  (3.10)  and  (3.11), 

S2  ■ p/ u , 

where  u is  given  by  (3.12). 


(5.5) 
A-  1,2,3 


(5.6) 


b . Reality  of  the  velocities  of  propagation  of  a wave 
Equations  (5.5),  together  with  the  equation 


l Va  " 0 (6.1) 

A 

are  equivalent  to  the  equations  (3.3)  and  (2.15)  of  [2]  which 

p 

were  used  to  calculate  the  possible  values  of  p/S  . It  was 
found  in  [2]  that  there  are  two  such  values  which  are  given  by 
the  roots  of  the  quadratic  equation 

(£)* •"(£)*“■ 0 ’•  (6-2) 

where 

a - 4(A^  + ...){(X^K2K3*  ...) 

♦ \ - *l)\*  ...J} 

♦ - X^)2M1+  ...] 

♦ 16n2n2n2(X*  - \2)2(\2  - X^)2(X^  - x2)2(WxlW22  - W^)  , 

p O p p ^ 

6 « ZU^CXJn^  X^np  ♦ ...] 

♦ [(A2-  \23)2n22n23 M1+  ...]}  , 

and  Ka,  Ma  (A  ■ 1,2,3)  are  defined  by 


V 


AAW2 


M. 


2(K11* 


2AAK12* 


x!*hl«) 

A 22 ’ 


(6.4) 


In  (6.3),  the  dots  indicate  terms  obtained  from  those  given 
by  cyclic  permutation  of  the  subscripts  1,2,3  on  the  X's, 
n' s , K ' s and  M's. 

As  was  seen  in  [2] , the  necessary  and  sufficient  conditions 
that  the  velocities  of  propagation,  parallel  to  the  unit  vector 
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n , of  sinusoidal  waves  of  infinitesimal  amplitude,  be  real  are 
ct  > 0 , 0 > 0 , B2  - 4a  > 0 . (6.5) 

In  [2]  the  following  expression  was  also  given  for  02  - 4a  : 

B2-  4a  « 4W2{[xJn!;(X2-  X2)2-  ...] 

♦ 2 [X2X2n2n2(X2-  X2)(X2-  X2)  + ...]} 

♦ 8W2{[(X2-  x2)2n2n2Mi{x2(x2.  X*)  (nJ  + n|) 

♦ Xl(Xl’  ^2)(ni+  n3)}]  + •••} 

+ 4{[(X2-  X2)2n2n2 M1+  ...]2 

♦ 16n2n2n2(X2-  X|)2(X|  - X2)2(X2-  X2)2(W22-  W^W^)  } . 

(6.6) 

By  substituting  in  (6.6)  the  expression  for  n given  in  (4.1) 
and  the  expressions  (6.4)  for  KA  and  MA  , it  can  be  shown, 
after  lengthy  algebraic  manipulation,  that 

02  - 4a  - 16(P2+Q2)  , (6.7) 

where  P and  Q are  defined  by  (4.4).  „ with  (4.5).  Thus, 
the  relation  02  - 4a  > 0 in  (6.5)  is  automatically  satisfied. 

In  a similar  manner  it  can  be  shown  that 

6 - 4R  . (6.8) 

The  validity  of  (6.7)  and  (6.8)  can,  of  course,  be  verified 
pragmatically  by  choosing  arbitrary  numerical  values  for  W , 

W2*  wn»  wi2»  W22  and  arbitrary  values  for  aA,  bA,  nA,  XA 
which  satisfy  the  relations  (4.1)  and  (3.4).  Then  B2  - 4a 
and  0 can  be  calculated  by  using  (6.6)  and  (6.8)  and  P2+Q2 


Id. 


and  R can  be  calculated  from  (4.4)  and  (4.5). 

From  (t>.7)  and  (6.8)  it  follows  that 

a - 4(R2-P2-Q2)  . (6.9) 

Equations  (6.7),  (6.8)  and  (6.9)  imply  that  the  conditions  (6.5) 
for  reality  of  the  wave  velocities  are  precisely  the  same  as  the 
conditions  (4.8)  for  the  incremental  shear  modulus  to  be  positive. 

It  was  shown  in  [3]  that  if  we  restrict  the  direction  of 
propagation  of  the  wave  to  lie  in  a principal  plane  of  the  under- 
lying pure  homogeneous  deformation,  then  the  necessary  and  suffi- 
cient conditions  that  the  velocities  of  propagation  be  real  are 

KA  > 0 and  flx  - * r;)MA+  Ka  > 0 (A -1,2, 3).  (6.10) 

It  was  shown  in  [Ij  that  these  are  also  the  conditions  that  the 
incremental  shear  modulus  be  positive  if  the  plane  of  shear  is 
a principal  plane.  It  was  further  shown  in  [2)  that  the  condi- 
tions (6.10)  imply  that  the  relation  (6.5),,,  and  hence  the 
relation  (4.8)  , is  satisfied  for  arbitrary  choice  of  the  unit 
vector  n . However,  it  was  shown  in  [4]  that  they  do  not 
imply  that  the  relation  (6.5)^  and  equivalently  the  relation 
(4.8)2,  is  satisfied  for  arbitrary  n . 
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7 . Directions  of  polarization  of  the  wave 

From  (6.2)  it  follows,  with  (6.7)  and  (6.8),  that 

~ " y{6±  (B2-4a)’5}  - 2 { R ± (P2*Q2)’5}  (7.1) 

S2  1 

for  a wave  propagating  in  the  direction  of  n . We  see  from 
(5.6)  that  p/S*"  is  the  value  of  the  incremental  shear  modulus 
for  shear  in  the  direction  of  polarization  of  the  wave,  the 
plane  of  shear  being  parallel  to  n . 

From  (4.10)  it  is  seen  that  the  expressions  on  the  right 
hand  side  of  (7.1)  are  the  stationary  values  of  the  incremental 
shear  modulus  for  specified  n . It  follows  that  the  directions 
of  polarization  of  the  sinusoidal  waves,  which  can  be  propagated 
parallel  to  n , are  the  directions  for  which  the  incremental 
shear  modulus,  for  directions  of  shear  perpendicular  to  n and 
•plane  of  shear  parallel  to  n , has  stationary  values. 
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